home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2001 May / SGI IRIX Base Documentation 2001 May.iso / usr / share / catman / p_man / cat3 / complib / clarfg.z / clarfg
Encoding:
Text File  |  1998-10-30  |  1.9 KB  |  67 lines
  1.  
  2.  
  3.  
  4. CCCCLLLLAAAARRRRFFFFGGGG((((3333FFFF))))                                                          CCCCLLLLAAAARRRRFFFFGGGG((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CLARFG - generate a complex elementary reflector H of order n, such that
  10.      H' * ( alpha ) = ( beta ), H' * H = I
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE CLARFG( N, ALPHA, X, INCX, TAU )
  14.  
  15.          INTEGER        INCX, N
  16.  
  17.          COMPLEX        ALPHA, TAU
  18.  
  19.          COMPLEX        X( * )
  20.  
  21. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  22.      CLARFG generates a complex elementary reflector H of order n, such that
  23.                 (   x   )   (   0  )
  24.  
  25.      where alpha and beta are scalars, with beta real, and x is an (n-1)-
  26.      element complex vector. H is represented in the form
  27.  
  28.            H = I - tau * ( 1 ) * ( 1 v' ) ,
  29.                          ( v )
  30.  
  31.      where tau is a complex scalar and v is a complex (n-1)-element vector.
  32.      Note that H is not hermitian.
  33.  
  34.      If the elements of x are all zero and alpha is real, then tau = 0 and H
  35.      is taken to be the unit matrix.
  36.  
  37.      Otherwise  1 <= real(tau) <= 2  and  abs(tau-1) <= 1 .
  38.  
  39.  
  40. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  41.      N       (input) INTEGER
  42.              The order of the elementary reflector.
  43.  
  44.      ALPHA   (input/output) COMPLEX
  45.              On entry, the value alpha.  On exit, it is overwritten with the
  46.              value beta.
  47.  
  48.      X       (input/output) COMPLEX array, dimension
  49.              (1+(N-2)*abs(INCX)) On entry, the vector x.  On exit, it is
  50.              overwritten with the vector v.
  51.  
  52.      INCX    (input) INTEGER
  53.              The increment between elements of X. INCX > 0.
  54.  
  55.      TAU     (output) COMPLEX
  56.              The value tau.
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.